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Supermassive Black Holes as Giant Bose-Einstein Condensates 
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PACS 04 . 70 . Bw - Classical black holes 

PACS 04 . 20 . Cv - Fundamental problems and general formalism 

PACS 04 . 20 . Jb - Exact solutions 

Abstract. - The Schwarzschild metric has a divergent energy density at the horizon, which 
motivates a new approach to black holes. If matter is spread uniformly throughout the interior of 
a supermassive black hole, with mass M ^ M* = 2.34 W^Mq, it may arise from a Bose-Einstein 
condensate of densely packed H-atoms. Within the Relativistic Theory of Gravitation with a 
positive cosmological constant, a bosonic quantum field is coupled to the curvature scalar. In the 
Bose-Einstein condensed groundstate an exact, selfconsistent solution for the metric is presented. 
It is regular with a specific shape at the origin. The redshift at the horizon is finite but large, 
z ~ 10^*M*/M. The binding energy remains as an additional parameter to characterize the BH; 
alternatively, the mass observed at infinity can be any fraction of the rest mass of its constituents. 



On the basis of the Schwarzschild, Kerr and Kcrr- 
Newman metrics, it generally believed that black holes 
(BHs) are singular objects with all matter localized in 
the center or, if rotating, on an infinitely thin ring. Re- 
cent approaches challenge this unintuitive assumption and 
consider matter just spread throughout the interior. [1~3] 
Here we shall follow this line of research. To start, let us 
just look at some orders of magnitude. For solar mass neu- 
tron stars it is known that the density is about the nuclear 
density. Solar mass black holes are about 10 times smaller, 
thus 1000 times more compact. Clearly, this begs for a 
Quantum Chromodynamics description in curved space. 

We shall focus on the other extreme, superheavy BHs. 
They occur in the center of each galaxy and weigh about 
-^^BH = 0.0012 Mbuigo- [4] Let us assume that they consist 
of hydrogen atoms and that mass and particle number are 
related as M = vNtu-h with some v < 1. If we neglect 
rotation, the radius is i? = GMjc? (see below). We may 
compare the BH density ZN /AttR^ with the one of densely 
packed, non-overlapping H-atoms, that is, with the Bohr 
density tt-b = 3/47rao, with oq = 0.529 A the Bohr radius. 
This yields a mass M* = c^{aQ/G^mY{Y/'^ (we take v = 1 
here), which lies in the range of observed supermassive 
black holes, AU ^ 4.66 lO^^kg = 2.34 10^ Mq. 

Next comes the question how matter can withstand 
the enormous pressure normally associated with such high 
densities. It was proposed originally by Sacharov that the 
vacuum equation of state p = —p could describe matter 
at superhigh densities. [5] Laughlin and coworkers assume 



that matter near the horizon could be in its Bose-Einstein 
condensed (BEC) phase, modeled by the vacuum equa- 
tion of state. [1] Dymnikova considers BHs obeying it in 
the interior, which, however, have one or two horizons. [2] 
Mazur and Mottola take the BEC idea over to the interior, 
and investigate a "gravo-star" , of which the interior obeys 
the vacuum equation of state, and which is surrounded by 
a thin shell of normal matter having the stiff equation of 
state p = +p. This solution is regular everywhere. [3] 

We shall demonstrate that a supermassive BH can ex- 
ist as a self-gravitating hydrogen cloud, in a Bose-Einstein 
condensed phase. We study the problem in a series of im- 
proved starting points: assume a stiff equation of state, 
self-consistently solve a quantum field coupled to the cur- 
vature scalar, first for a uniform groundstate wavefunc- 
tion and next for a space-dependent one. Hereto wc have 
to employ the Relativistic Theory of Gravitation (RTG), 
which reproduces all weak gravitational effects in the solar 
system [6,7] as well as the ACDM cosmology [8]. 

We consider a static metric with spherical symmetry, 



ds^ = Uir^dt^ - V{r)dr^ - W^{r){de^ + sin^M^^). (1) 

The gravitational energy density arises from the Landau- 
Lifshitz pseudo-tensor, [9] generalized to become a tensor 
in Minkowski space. [8, 10] For (1) it takes the form 
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Let us start with the General Theory of Relativity (GTR). 
The Schwarzschild metric reads in the harmonic gauge 



1 _ 2M _ r - M 



Ws = r + M. 



(3) 



(We put G = c=?i=l.) It is singular at the horizon 
rfi = M and involves the gravitational energy density 
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Its quadratic divergence at r^ presents a hitherto over- 
looked peculiarity, that induces a negative infinite contri- 
bution to the total energy. For this reason, we shall switch 
to RTG with matter not located at the singularity r = 0, 
but just spread out within the horizon. 

Quantum field theory of Bose-Einstein condensed black 
holes. Let our H-atoms be described by a bosonic field 



^('''O = X!^*'^»w 
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energy density Eq. (2). [8] It extends the Hilbert-Einstein 
action with the cosmological term and a bimetric coupling 
between the Minkowski (7) and Riemann (g) metrics. 
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(For pbi = it is just a field theory for GTR.) One has 
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with i = r,9, (j). The value pbi = PA is imposed to have a 
Minkowski metric in the absence of matter. One may fix 
them to the observed positive cosmological constant [8]. 
However, historically the opposite choice pA < was con- 
sidered and the cosmological data were described by an 
additional inflaton field. [7] So the sign of pbi is not known 
yet; We show that solving a realistic black hole settles this 
issue. The new point of RTG is that 500 = U can be very 
small. Despite the smallness of pbi, the phi/U terms be- 
come relevant near the horizon [6,8], bringing 



where i = {n, £, m}, [ai, dA = Sij and eigenfunctions factor 
as tpi{r) = 4'n(j')Yim{0, <p)- The rotating wave approxima- 
tion then leads to the Lagrangian [11] 



L^.t = g'^'d^^^dj - (m2 + ^R)^^ - -^^ V^ (5) 

For a field in curved space the renormalization group gen- 
erates a coupling to the Ricci curvature scalar R. [12] 
Its strength ^ is for now a phenomenological parame- 
ter. The dimensionless coupling A — 8nr^cg/h with 
g = Awh as/m models the two particle interaction by the 
scattering length. For hydrogen in flat space one has [13] 



as = 0.32 flo singlet state, a^ — 1.34ao triplet. 

We shall continue with the singlet value A = 0.81 10^. 

With *o = (2£;o^o)^/^V'oe"*-^''* for Nq groundstate 
atoms, the relativistic Gross-Pitaevskii equation reads 



^ 5^V=55^''9,+m2+ei?+ 4^1*0 = 0. (6) 



9 ■ " ' ^Eo 

A homogeneous ground state, 'i'o{r,t) = ^Poe^*^"* 



curs when 



-f +m^+ei?+^l*o|^ = 0. 



(7) 



We focus on RTG, which describes gravitation as a field in 
Minkowski space [6,7] and posesses the same gravitational 
energy momentum tensor and thus also the gravitational 



R = -87rTtot - M-P + Pr+Pe+P^ + j^)- (10) 
To start, let us consider the stiff equation of state 

1), p.=p^]-p.{^-l), (11) 
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where pc = 3/327rAf^. For Uc = it is the vacuum equa- 
tion of state p = —p =const. One gets 



R = 8tt 
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Eq. (7) has a solution due to the ^i?-term. The constant 
and 1/U terms term imply. 



e = eo(i 
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:), E^ = 87raPcUc + Phi), (13) 



respectively. The dimensionless parameter 



Co = -m^M^ 



1.8010'^N ^ 



(14) 



appears to be large, but since R ^ 1/M^ the combination 
^R is just of order ?Ti^, making its effect of order unity. 
The metric can be solved as below; we omit details. 

Self- consistent field theory. Rather than imposing an 
equation of state, the material energy-momentum tensor 
T^'' should be derived from first principles, i.e. from the 
quantum field theory for the H-atoms. Its energy density 
reads, if we exclude the effect of the ^i?-term. 
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The pressures {p^^ , p]p , pT) have this shape with signature For the Schwarzschild black hole the horizon occurs when 

(+H ), (H \ ), and (H 1 ), respectively. A4 — M for W = 2M. Concerning the outside metric, we 

Spherical symmetry will imply that p™ — p^ = p™. For will be close to that situation. This implies again that a 

a uniform groundstate Pm is isotropic, mass M corresponds to ^ = 2m?M^ /3. 



. Let us introduce the 'Riemann' variables x and y by 

{Pm,Pm) = ^{^ ± ^)\^l\ ± ^. (16) 



a: = — , y=Vl-a;2, (26) 

They consist of a vacuum part p = —p = const, and a stiff 

part p = +p ^ 1/U, the types studied in [3] and above, so that U = 2p,'^x^. The pbi terms in (8) violate general 

In the non-relativistic {Eq = m) and flat space {U = 1) coordinate invariance and impose the harmonic gauge, 

limit, they reduce for A = to pm — mc^ |\I/g| and pm — 0- 

Because of the fi?-term in (5), the Einstein equations 

. U V W ArV 
embody a direct backreaction of matter on curvature, 1-4 — (27) 

G^'^ ==:87r(r;^'^ + r;;"' + T^i'')-167rC(V'^)G^''. To connect U V W W^ ' 

to the standard notation, G^'' = SnT^. we define T^^- by ^-^^^ (22), (23) and M = \M{x + x^) from (25), it brings 

l + B 
with direct backreaction strength of matter on the metric 



TZ^ " , ^„ ""' =T^" + r^" + T,^7, (17) 2x' 2x" 2xx' Ax' Srx 

l + B ' — 



X x' 1 — X^ X x^{l — x^) 

Going to the inverse function r{x) makes it linear. 



B = 16<(VV) = ^4!^- (18) x'{l-x'y' + x{3-4x'y=Ar. (28) 

Ea 

T^ \ r, ii i 1 r 11 r /i ^\ Thc solution is thcn remarkably siuiplc, 

_bor A = the curvature scalar follows from (16) as .; j- i 

^^ 87r / E„\^l\+p^, ^"^'l-PglV (19) r^n{l + ^)x^''\l + y)'^\ (29) 

1 + -B \ U Eq J 

Solving Eq. (7), we find two relations and a consequence, (The second independent solution with V5 ^ -V5 is sin- 

gular.) This determines the metric function V, 

5 = 1, So' = 8<pbi, l*ol = ^ = ^- (20) 2VK'2 5M2 

8< ^0 V = l^ = ^x4-2^(l + y)2^. (30) 

(The GTR situation, reached by taking pbi -^ first, 

would not allow a meaningful solution.) The first iden- Putting these results together, it now follows that 

tity expresses a 100% direct backreaction of matter on the 

metric. This motivates to introduce the parameters [81 „ _ „ _ /'qi \ 

fi = ^/l6^^pbl = \/2A, p = pM = 7.90 10"^^^- (21) ^s the 1/U terms cancel due to the relation E^ = Sir^pbi- 

So, after all, we reproduce the vacuum equation of state. 

Instead of searching a finite U, as for boson stars, [14] we ^o normahze [foj, we need the 3d volume element in the 

assume a very smaU U with C/(0) = 0, coded by v, future time direction, dE'^ = drd0d(|)n^' y/^^ = S^dV, 

j^ set by the timelike unit vector n^ — 6q/\/U and g^ = 

U = 2^^^"^^^ *^^^'' -VW^ sin^ 9. This results in dV = drdn^/V/UW'^. 

The general inner product [12] (i/'i,'02) = — ■j/dS'^x 
{'>pidfj,'>p2 ~ '02^A''0i) defines the orthonormality 



In terms of the mass function A4 (r) , defined by 



(V'» , ^j ) = [E, +E,) dV ^.^* = 5^, . (32) 
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the 00 and 11 Einstein equations take the form 

W^2MU' M , , ,,,, With dV = d2/dr!8M3/(M^') it yields 

2UW^ W W^ 

The Ansatz (22) solves them and yields, due to (20), 1^21 _ nrp ly l/,2| _ ^o'^P' /oo\ 
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In Eq. (13) the correction term is of order 
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(34) 
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These corrections seem relevant for BH's with masses 
M ~ 2.75 lO-^M* ~ 645 M©. Much less below M* the hy- 
drogen atoms will get ionized, calling for fermionic fields 
for protons and electrons, which by a BCS pairing can 
again undergo a BEC transition. This BCS-BEC scenario 
is beyond the aim of the present paper. 

The exterior. At the horizon r/j w ri , y/j <C 1 one has 
5 „, ... „„ 1/^ 



U 



>yh- (35) 



2/2^, V^'^, W^2M, W 2 

We have to connect this to the vacuum solution outside 
the BH. Well away from matter, the harmonic constraint 
brings the Schwarzschild shape (3), where M = A4{rh) is 
the mass, essentially as observed at infinity. The values 
(35) are far from Schwarzschild's, even when r is near 
M (e.g., Wg = 1). The problem nevertheless appears to 
be consistent. Near rh we need the deformation of the 
Schwarzschild metric which regularizes its singularity due 
to the bimetric coupling. [6, 7] An elegant scaling form 
for small p, was presented by us, [8], 
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W = 



U = rye?, (36) 
2M 



??(l + e«) 



{\2' 



1 — T]ei - yU^(C + Wo) 



Here ^ is the running variable and 77 a small scale. For 
Tjc^ = 0(1) it coincides with the Schwarzschild solution. 
Matching with the interior appears to be possible. 
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Vbfl, r; = — M, M/' = e«^ + ^ ^ ^V5m, (37) 
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implying jjh = 3/1. Taken together, the three regimes, 
interior, horizon and exterior, provide an exact solution of 
the problem. At the origin it exhibits the singularities 
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where 7^ = ^(VS + 1) is the golden mean. But if we take 
W as the coordinate, we have in the interior the shape 



ds2 ^ -^i^WMt^ 
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which is regular at its origin, with the term 2dW^ coding 
the above singularities. 

We may rewrite the exterior solution by eliminating ^, 
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This describes the free space region r > M, where 2p.^ < 
U < I + 0(fi). At the cosmic scale r ^ l//i Newton's 
law picks up the Yukawa-type factor cos iir, due to the 
tachyonic nature of gravitation in RTG with pbi > 0. [8] 

The interior shape can also be expressed with U as run- 
ning variable, where it lies in the range (0,2/22). j^^^ ^^ 
Eq. (29) it also holds that 



W = 2Mx 
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'2U, 




(41) 



The locus and the metric function V are given by (29), 
with ri « M, and (30), respectively. 

Contributions to the energy. With the weight dV 
given below (32), the standard expression for the energy, 
/ dV/O scales as 1/p, and even diverges logarithmically at 
r = 0. However, in RTG the energy is determined by Eq. 
(2). At the origin it diverges as r^^"^, which is integrable. 
The gravitational energy inside the BH reads 



t.^grav, int — ^^ 



drr^t°° = -842.898 M. 



The total energy density reads 6°° = t°° + VW^ptot/r'^- 
We can calculate the material and the bimetric energy. 
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169.431 M, 



= 686.466 M. 



Together they make up for 
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^grav, int 

The energy density in the skin layer first has a large posi- 
tive and then a large negative part, due to the term r^V' . 
The integrated effect is obtained easily since, in the formu- 
lation of the Einstein equations in Minkowski space, the 
total energy density is a total derivative, [8] The region 
r > M thus yields t4xterior = — 12M. Together with the 
interior it confirms the total energy U — Mc^, expected 
from the decay of the metric, goo = 1 ^ 2GM/c^r. 

Non-uniform groundstate. Till now we assumed that 
the groundstate wavefunction has a constant amplitude, 
and drops to zero at the horizon. Clearly, this cannot be 
exact. Taking ^0 -^ ^o(^), we first take into account that 
in deriving the Einstein equations, partial integrations are 
to be performed. This brings derivatives of i? ~ l^oP) 
and induces an extra term T^'^, 
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Eq. (43) now leads to a total energy momentum tensor Expressing the shapes (50) in y, we have 
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The Gross-Pitaevskii equation (6) reads in terms of B 
B" 
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We can now first verify that the total energy momentum 
tensor is conserved due to the harmonic condition (27). 
The singular term B /U also drops out from (47) for 
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This was already used in (46) to cancel the 1/U terms. 
With 5 = 2m'^M'l/?, and /z = fiMi it implies again _Eo = 
Jlm/^/ji. For very small B <^ 1/S, ^ 10~^^ there will be an 
exponential fall off, B ^ exp(— \/5?Tir), so the horizon is a 
fraction of the Compton length thick. In this narrow range 
p, Pr and p_L vanish smoothly. In the regime B » 1/^, 
Eq. (47) simplifies and actually reduces to Eq. (7), 
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Here we can consider B as vanishing sharply, B ^ r^ ^ r, 
so that Tg'^ y^ 0, keeping the ultimate exponential tail and 
decay of p and the Pi in mind. Eq. (46) then brings 
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Let us first return to the interior where U — 2p,'^x^, V - 
8Mfx'^/y^ and W = 2Mix. Eq. (49) can be written as 
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yBy + (1 - X)B = 1. 
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To understand the structure of the problem, we again take 
A = 0. Then for any A there is the solution 



B{x) = l + Ay = l + A\J\-x'^, 
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{2B + yBy+4), 
{2B-yBy). (53) 



Surprisingly, their A-dependence factors out, keeping a 
vacuum equation of state p = —p = S/GAnAlf, so (52) is a 
non-uniform, exact solution of the same metric. The hori- 
zon _B = is now located at r^ > ri where yu = —1/A. 
[Eq. (29) continues to negative y « V5{ri — r)/Ari for 
r > ri]. However, a problem shows up with the match- 
ing, since W'{rh) ^ ~1/^ cannot be of order jl ~ 10~^^ 
anymore. We thus have to deviate from the exact solu- 
tion, which leads in general to a numerical problem. An- 
alytically, this question can be considered for large A, by 
adding \/A^ corrections to previous solution. Expanding 
in 1/A at fixed s = A\fh{r jr^ — l)/4, we arrive at 
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The additional terms, found to be 



i>\ = ho + &11S, 

Vl = 1^10, 



ui = uio + 4wii ln(2 — s), 

wi = wiQ — wii ln(2 — s), (55) 



produce an anisotropy, p ^ —pr ^ —pi_ 7^ P- The horizon 
S = is located at s^ = 1 + ^/5/2A + 0(1/^2), where 
W' = V^{wn - 1)/2A. Clearly, W' -^ p from (37) can 
be attained by tuning wu = 1-1- 3Ap, + 0{1/A). The 
maximum of VF at s = in the absence of the wi term 
has now been shifted to the horizon, which shows that the 
problem has a proper solution, with A remaining a free 
parameter. U has a maximum at s = 1 — \/3. The mass 
seen at infinity, M = M{rh) — ^W{rii), coincides with 
Ml to order 1/^2. At the horizon we can fix ri from (36), 



n = Ml 



1- 



-±j + i|(2 + ,n2,') + 0(-l: 



(56) 



Properties of the solution. The groundstate occupa- 
tion number becomes upon neglecting the l/yl2 corrections 



/M f^ 
dV|*o|'-2V3— / dy(l 



+ Ay). (57) 



We may write the two leading orders as 
M = i^Nom, 



V3{2 + A) V3A' 
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Clearly, the energy Mc^ of the BH can be any fraction of 
the rest energy Nomc^ of the constituent hydrogen atoms 
If V starts at a value z^c < 1, our BH is likely approached 
in an explosive manner, possibly related to jets of quasars. 
We found the sharpness of the horizon to be a fraction 
of the Compton length of H. On a much larger scale i?grav 
there is near- horizon growth of the metric functions U/t] « 
rjV ss e^, taking place in the millimeter range, 



dr 4 

^grav - 37 = 2?7M =—=flM = 4., 



10-3 _^. (59) 



It is a realistic value, small compared to the size of the 
BH, and still large compared to the Bohr radius. 

The interacting situation. If the nonlinearity A is rele- 
vant, a numerical solution is called for. Assuming the same 
leading order behaviors near r = 0, the above structure 
survives. Mass and particle number remain independent 
parameters. This analysis remains as a task for future. 

Conclusion. We have questioned the general wisdom 
that static BHs have all their mass in the center and can- 
not be described by present theories. Numerical estimates 
show that a picture of closely packed H atoms naturally 
applies to the supermassive BH's in the center of galaxies, 
M ^ Mi, ^ 2.34 lO^M©. We present within the Rela- 
tivistic Theory of Gravitation (RTG; a colloquial term is: 
Not-so- General Relativity), an exact solution for a BH, of 
which the interior is governed by quantum matter in its 
Bose-Einstein condensed phase. Its density decays alge- 
braically in the bulk and exponentially near the horizon. 
This solution is matched with the Schwarzschild metric, 
which near the horizon is deformed in RTG. Powerlaw 
singularities occur at the origin, that get absorbed in the 
Riemann description of the metric. Elsewhere, the solu- 
tion is regular. The redshift at the horizon is finite, though 
of the order 1/jl ^ 10^^M*/M. To specify a BH requires 
not only the mass (and, in general, charge and spin), but 
also the rest energy of the constituent matter. 

Our BH is a quantum fluid confined by its own gravi- 
tation. In the interior, time keeps its standard role. No 
Planckian physics is involved; Hawking radiation is absent 
and Bekenstein-Hawking entropy plays no role. 

Our BH has one "hair" . As one would expect for a 
classical theory of gravitation, when the quantum matter 
in the BH has reached a certain groundstate, the classical 
metric allows the system still to go to a lower energy state. 
Indeed, the passage of celestial bodies will induce oscilla- 
tions in the metric and emission of gravitational waves, 
which, upon re-equilibration, increase the binding energy, 
finally up to 100% of the rest energy of its constituents, 
Nmc^. This property may explain the enormous jets and 
energy output of quasars and also be responsible for very 
high energy cosmic rays, E > A lO^^eV. [15] 

An important question is whether formation of realistic 
supermassive BHs brings the matter indeed in or near the 
Bosc-Einstcin condensed groundstate. Also the stability 
of the solution needs to be studied. It also remains to 



be seen whether the phenomcnological value for the pa- 
rameter ^ has a microscopic underpinning. Extension to 
finite temperatures, not presented here, will exhibit a T"^' ^ 
fraction of thermal atoms. 

Calculation of the normal mode spectrum may lead to 
predictions that deviate from the ones of GTR; this spec- 
trum may be observed in the foreseeable future. 

If we apply Eq. (48) to Logunov's case pbi < 0, it follows 
that ^ < 0, so i? < due to (18). To avoid a singularity in 
e.g. (46), a lower bound -B(O) > —1 will be required. As 
already indicated by the exact solution (52), B will then go 
to zero at some point well below M. This prevents fitting 
to the external metric and excludes our BH solution. 

We failed to apply our approach to GTR, technically 
because it lacks compensation for the singular 1/U terms. 
If no other solution exists for the considered physical sit- 
uation, GTR must be abandoned and replaced by an- 
other theory, RTG being the first candidate. In view of 
its smaller symmetry group, this may have far reaching 
consequences for singularities in classical gravitation and 
for quantum approaches to gravitation, while Minkowski 
space-time needs no quantization. 



The author has benefited from discussion with Steve 
Carlip, Kostas Skenderis and Bahar Mehmani. 
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